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ABSTRACT 

Mergers are predicted in all cosmological models involving Cold Dark Matter to 
be one of the dominant channels whereby galaxies accrete mass. In this paper we 
present a detailed analysis of predicted galaxy-galaxy merger fractions and rates in 
the Millennium simulation and compare these with the most up to date observations 
of the same quantities up to z ~ 3. Wc carry out our analysis by considering the 
predicted merger history in the Millennium simulation within a given time interval, as 
a function of stellar mass. This method, as opposed to pair fraction counts, considers 
mergers that have already taken place, and allows a more direct comparison with the 
observed rates and fractions measured with the concentration-asymmetry-clumpiness 
(CAS) method. Wc examine the evolution of the predicted merger fraction and rate 
in the Millennium simulation for galaxies with stellar masses ~ 10^ — 10^^ Mq . 
We find that the predicted merger rates and fractions match the observations well for 
galaxies with Mi, > 10^^ Mq at z < 2, while significant discrepancies occur at lower 
stellar masses, and at z > 2 for Mi, > 10^^ Mq systems. At z > 2 the simulations 
underpredict the observed merger fractions by a factor of 4-10. The shape of the pre- 
dicted merger fraction and rate evolutions are similar to the observations up to z ~ 2, 
and peak at 1 < z < 2 in almost all mass bins. The exception is the merger rate of 
galaxies with > 10^^ Mq , which remains high at z < 1.5. We discuss possible 
reasons for these discrepancies, and compare different realisations of the Millennium 
simulation to understand the effect of varying the physical implementation of feed- 
back. We conclude that the comparison is potentially affected by a number of issues, 
including uncertainties in interpreting the observations and simulations in terms of 
the assumed merger mass ratios and merger time-scales. The differences between the 
observations and simulation results might also be due to problems in the modelling 
of star formation in the simulation, which produces redder and less biased galaxies 
than observed, particularly for galaxies with stellar masses Mi, < 10^^ Mq Finally, 
our findings may also be related to other CDM problems, including the the lack of 
massive galaxies with Mi, > 10^^ Mq at z > 1, and a lack of merger events between 
lower mass galaxies. 
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1 INTRODUCTION 

Theoretically, the dominant and largely successful paradigm 
for understanding and modelling the universe is based on 
the theory of Cold Dark Matter (CDM). In this scenario, 
the first objects that form in the universe are low mass 
haloes that subsequently merge togeth er to form more mas - 
sive structures as time progresses (e.e:.. IWhite fc Reeslll978l ). 
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The hierarchical mass assembly by merging is in fact the 
cornerstone of CDM-based simulations and can be empiri- 
cally tested, for example, by investigating the rol e of mergers 



in the formation and evolution of galaxies (e.g. KauffmannI 



1 19931 : iBerrier et al1l2006l : IWetzel. Cohn fc Whitell2008D . 

Galaxy formation and evolution is certainly driven in 
part by galaxy mergers. There is no doubt that in the lo- 
cal universe ther e are examples of gala xies merging with 
each other (e.g.. [Pe Propris et al.|[2007^ ■ which eventually 
will evolve into a system more massive than either of the 
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progenitors. During this process the physical evolution of the 
resulting remnant galaxy will also change, as mergers are po- 
tentially the dominant method for triggering star formation 
and producing feedback by ejecting, or heating, gas (e.g., 
ICox et al-llioM) . Mergers can also spur on the formation of 
black holes and trigger AGN a ctivity, as both observations 
and theory have shown (e.g., iHopkinj l2008l : iBundy et al.l 
|2008| ). While galaxy mergers are seen in the local and dis- 
tant universe, the exact role of mergers in the formation and 
evolution of galaxies over cosmic history is still uncertain. 

There are two primary methods for tracing the 
merger history of galaxies in observations: morphologi- 
cal identific ation techniques a nd the close galaxy pair 
method (e.g. iPatton et aLll2000t ). Examples of morphologi- 
cal techniques are the concen tration-asymmet ry-clumpiness 
method (CAS, hereafter, e.g. 'Conselice 2003a') and the Gini- 
M20 method (Lotz, Primack & Madau 2004). The CAS 
method identifies galaxies undergoing mergers o n the ba- 
sis of observed morphological properties (e.g., IConselic3 
l2003al ) and is a non-parametric technique for measuring the 
shap es and structures of galaxies on resolved CCD images 
(e.g.JConselice. Bershadv fc Jangren.2p00a: Bershadv et al.l 
I2OOOI : IConselice et al.ll2002l : IConseliciboOSal ). AU method- 
ologies have been demonstrated to be internally consis- 
tent with each other jPe Propris et al ]|2007l :[ Conselice et al.l 
l2008l : iKitzbichlcr fc Whitd 120081 '). However, since each 
method identifies galaxies at a different stage during the 
merger, overall the galaxy populations investigated by each 
method can be substantially different. The main uncer- 
tainty in all methods resides in the time- scale Tm used to 
identify galaxies u ndergoing mergers (e.g.. IConselicell2006al : 
iLotz et al]l2008bl ). Currently, the most likely time-scale to 
which the CAS and pair methods are sensitive is believed to 
be Tm ~ 0.4 - 1 Gyr. 

In the last decade a r elatively large sample of galaxy 



mergers has been collated (IPatton et al. 
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De Propris et all 12007: Lot_z_ct_,al. 2008a,; Conselice et all 



20081 : iPatton fc Atfieldl l2008i : iBluck et al.l I2OO9I ) and the 



merger history of galaxies has been traced with sufficient 
accuracy to perform a detailed comparison with theoretical 
models. Mergers are just now beginning to be identified at 
z > 2 through either close galaxy pairs (e.g., Bluck et al.l 
20091'). or structural methods (e.g., Conselice et al.ll2003bl : 



Conselice et~alll2008l : iLotz et al.l l 2008bl '). What is becoming 



clear is that the merger fraction increases at larger look back 
times (higher redshifts) . This increase can be parameterised 
as eit her a power-law, or a cornbined power-law/exponentia l 
(e.g.. lCarlbe"r3ll99Gl : IConselicdl2006 al: IConselice et al ] |2008h . 
These measured merger fractions in the early universe can 
be converted into merger rates up to z ~ 3, and based 
on these we can attempt to determine the role of major 
mergers in the formation of galaxies. 

In this work, we retrieve the predicted galaxy merger 
history e volution from t he M illennium-based semi-analytic 
model of lBertone et al.l |2003) as a function of stellar mass. 
We then compare the model predictions with a selection 
of observational results and examine separately the merger 
history of galaxies with stellar masses in the intervals 10^ 
M© < Kh < 10^° Mq , > 10^° Mq and > lO" Mq , 



at 2; < 3. Given the uncertainty in the observational results 
introduced by the time-scale Tm, in this paper we assume 
Tm ~ 0.4 Gyr as our fiducial value and we investigate how 
results can change for a longer time-scale of Tm = 1 Gyr. 

In general, we find that the predicted merger history 
varies as a function of stellar mass, redshift and merger mass 
ratio. The predicted major and minor merger rates, that 
is the number of mergers per unit volume and unit time, 
increase with lower stellar masses. We also find that the 
predicted merger rate increases with redshift and reaches a 
plateau at 2 < 1 for galaxies with M* > 10^^ Mq . The 
predicted merger rate for less massive galaxies decreases 
with redshift. At z < 1.5 we find a good agreement be- 
tween predicted and observed merger rates. However, the 
predictions for the merger fractions match the observations 
only for the most massive galaxies with Ali, > 10^^ Mq at 
z < 2. The predicted major merger fractions are about 10 
times smaller than the observational estimates for galaxies 
with M* > 10^° Mq at 2; ~ 0.5. We discuss the various 
ways in which the merger fraction predictions could be in- 
correct and conclude that the simulations might not predict 
enough mergers between galaxies with Af* < 10^^ Mq . This 
could potentially account for other problems with matching 
observations to semi-analytical models, including the fact 
that there are more massive galaxies at high redshift than 
predicted and the sat ellite problem (e.g.. lMoore et al]|l999l : 
IConselice et al.ll2007h . 

This paper is organised as follows. In 321 we describe 
the observational data used to compare with the theoretical 
predictions. In particular, we briefiy discuss our criteria for 
selecting the observational sample and some basic proper- 
ties of the galaxies in the sample, such as their morphology 
and stella r masses. jS] de s cribes the Millennium galaxy cat- 
alogue of lBertone et al.1 l|2007l ) and how the merger rates 
and fractions have been extracted from the simulation. In 
331 we present the model results and compare them to the 
observational data. Finally, fJSl discusses our findings and JHl 
summarises our results. 



2 MERGER HISTORY DATA 

In this Section we give a brief description of the data sets 
used for comparison with the Millennium simulation (Sub- 
sections 12.11 and I2.2|l and of how the merger fractions and 
rates are measured (Subsection 12. 3|) . 

We use results from a series of works, based on both 
the structural asymmetries of galaxies and pair counts, that 
investigate the merger his t ory fo r gala^xies as a function 
of stellar mass (IConselicel l2003al: iDe Propris et aP l2007l: 
Conselice et"al] l2008l : iBluck et al.1 l2009l : IConselice et all 
2009h . Interested readers should examine these papers for 
the many details involved. The observational quantities we 
present in this paper are measured using a standard ACDM 
cosmology with Ho = 70 km s~^ Mpc~^ and Q,m = 1 ~ fl\ 
= 0.3. This differs slightly from the cosmology used within 
the simulation, although the effects of this are minor. 

The first step to investigate the merger history of galax- 
ies requires selecting a robust sample of mergers. The ob- 
served merger history must be constructed using different 
data sets, as the ability to measure merger properties at dif- 
ferent redshifts requires different observing conditions and 
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telescopes. This is due, e.g., to the fact that to obtain a large 
sample of nearby galaxies to measure the local merger rates 
and fractions requires a wide area survey, while detecting 
galaxies at high redshift requires very deep exposures, typi- 
cally with the Hubble S pace Telescope (e.g. IConselice et al.l 
l2008l : lLotz et al.ll2008al ;). 



2.1 CAS-Selected Mergers 

At low redshifts, we use datasets that estimate the merger 
quantities through the CAS method. The basic idea be- 
hind structural methods is that galaxies have light distri- 
butions that reve al their past and present formation modes 
(|Conselicell2003al ) . The classical rest-frame optical CAS def- 
inition for determining whether a sys tem is undergoin g a 
merger is given by the two conditions (|Conselicel l2003al ) : 



A > 0.35 and A > S, 



(1) 



that is, the asymmetry A must be larger than a thresh- 
old value, and must exceed the value of the dumpi- 
ness S of the galaxy. These selection cr iteria have been 
tested against sa mples of nearby galaxies (IConselicel 2003a 



De Propris et al][2007i l and N-body simulations ( Conselicg 
2006al ). Tests at higher redsh i fts a re limited to either 



smaU s amples (IConsehce et al.l I2OO8I '). or to the GEMS 
survey (|jogee et all |2009| ). The CAS method empirically 
finds systems which are nearly all ma j or mergers (e.g. 



Conselice. Bershadv fc GaUaghed l2000bl: IConselicd l2003al : 



Hernandez- Toledo et al.l l2005l : IConselicg l2006al ) , although 



there are hints that in some cases the mergers identified 
by CAS co uld be contamina ted by mergers with a smaller 
mass ratio l|jogee et al.ll2009l ). 

For galax;ies at z ^ 0.05, which is effectively our 
nearby galaxy merger f raction, we use the res ults from 
iDe Propris et al.1 (120071 '). iDe Propris et al.1 (|2007l ') estimate 
the merger fraction using both the number of systems in 
pairs and their morphological structure as identified by 
CAS. At 0.2 < 2 < 1.2 we take the merger fractions mea- 
sured for a combined COSMOS and Extended Groth Strip 
sample bv IConselice et al] (|2009l ). who consider more than 
20, 000 galaxies with M* > 10^° M© . The merger quanti- 
ties of systems at 1 < z < 2 are from a combined sam- 
ple of Hubble Deep F ield (HDF) and Hub ble Ultra-Deep 
Field (UDF) galaxies (|Conselice et al.ll2008h . The compila- 
tion of observational results we use in this wo rk is in good 



agreement wi th other published results (e . g. Patton et al.l 
20021: iLin et a l. 2004; Kartaltepe et al.ll2007l : lLin et al.ll200a 



Patton &i Atfieldi2008 ). A detail ed comparison betwee n the 



different datasets is presented in IConselice et al] l|2009l ) and 
we refer the reader to that paper for more information. 

2.2 Mergers from Pair Methods 

At the highest redshifts we consider, we use estimates of the 
merger quantities based on pair counts. The pair method 
is a complementary and independ ent approach for find ing 
the merger history of galaxies (see IConselice et al.ll200^ for 
details) . 

In the pair count method, close pairs are defined as 
galaxy systems with a physical separation of less than 30 
h^o' kpc. Major merger pairs are further defined as systems 



with galaxies with magnitudes within ±1.5 of each other. 
The pair fraction, defined as the number of galaxy pairs di- 
vided by the number of galaxies within the sample, is not a 
direct measure of the merger fraction. Th e conversion of a 
pair f raction in to a merge r fraction (e.g. iDe Propris et all 
I2OO7I : IConselice et al.|[2008l ) is affected by large uncertain- 
ties and in particular suffers from uncertainties in the deter- 
mination of the time-scale for merging. We will discuss the 
time-scale to which the CAS structural method and the pair 
method are sensitive in the next Sections. This time-scale is 
also required to convert fractions to merger rates and to ef- 
fectively compare the observed results with models. In this 
work we derive merger fract ions and rates from pairs only at 
z>2 based on lBluck et"aLl lf200a) drawn from galaxy pairs 
identified in deep NIR imaging from the GOODS NICMOS 
Survey. 



2.3 Merger fractions and merger rates 

In this Subsection we define the main quantities we will be 
investigating in the following, namely the merger fraction 
and the merger rate, and we discuss possible biases in their 
estimation both in observations and simulations. 

The merger fraction is the number of mergers A'^m within 
a given redshift bin and stellar mass limit, divided by the 
number of galaxies A^t within the same redshift 2 and stellar 
mass bin: 



/n.(M. 



A^m(M. 



A"t(M*,2)' 



(2) 



The conversion of the merger fraction, that is the ratio of 
the number of mergers to the total number of galaxies, into 
a galaxy merger fraction /gm (|Conselicell2006al ) , which is the 
number of galaxies merging divided by the total number of 
galaxies, can be done through: 



2xU 
1 + U 



(3) 



If each pair produces a merger, the pair fraction is equivalent 
to the merger fraction, sans different time-scales, and the 
value of the galaxy merger fraction /gm is rou ghly twice 
that o f the merger fraction /m, defined in Eq. (|2| ( Conselicg 
l2006al ). 

The conversion of pair and merger fractions into merger 
rates requires knowledge of the time-scale Tm to which the 
CAS and the pair methods are sensitive. This time-scale is 
critically important, because it ultimately determines the 
values of the merger rates and fractions both in simula- 
tions and observations. Simulations in general, including 
the Millennium, have well-defined merger rates, that can 
be measured accurately. However, the merger fractions de- 
pend on the time-scale over which the number of mergers are 
counted. In fact, the fraction of galaxies that have merged 
within a given time interval is roughly proportional to the 
time-scale considered. However, the number of mergers per 
unit time is independent of the time-scale used. 

The merger rate, defined as the number of galaxies 
merging per unit time and per unit volume, which we denote 
as 5R(2) , is given by: 



K(2) = 



/gm(2) ■ ng,n(2) 



(4) 
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where rigm is the number density of galaxies within a given 
stellar mass range. Although the observational value of Tm 
is not well constrained and there are a few indications that 
it might not be the same for structural and pair methods 
iLotz et al.ll2008al : lKitzbichler fc White! l2008l ). in general it 
varies between about 0. 4 Gyr and 1 Gyr (e.g., IConselicg 
l2006al : iLotz et al.|[2008bl V In this work we assume two dif- 
ferent values for Tm, that is Tm = 0.4 Gyr and Tm = 1 Gyr, 
in order to take in to account at first order the possible sys- 
tematics involved with comparing data with models. The 
value of Tm = 0.4 Gyr is based on the sensitivity of the 
CAS method, which is believed to identify galaxies that have 
undergone mergers within approximately such a time-scale. 
The tim e-scale of 1 Gyr has been suggested by iLotz et ahl 
l|2008bl ) on the basis of results of hydro-dynamical simula- 
tions. 

Another issue to address while comparing models with 
simulations is the minor vs. major merger sensitivity of dif- 
ferent methods. Pair fractions require galaxies to have mag- 
nitudes within ±1.5 of each other, but do not provide exact 
information about the respective galaxy masses. Although 
it is usually assumed that the mass ratio in galaxy pairs is 
larger than about 1:4, there might be large variations from 
pair to pair. As discussed in Section 12.11 the CAS method 
identifies galaxies that have already merged. This implies 
that we do not know the progenitor masses, although sim- 
ulations predict that minor mergers with mass ratios 1:10 
or less do not produce signifi c ant structural a symmetries 
iHernandez-Toledo et a"l]|2005l : IConselicell2006ah . The lim- 
ited information available suggests that the CAS method 
only picks out major mergers with a mass ratio of about 1:4 
or greater, but does not e xclude contamina tion from mergers 
with smaller mass ratio l|Jogee et al.ll2009h . 



3 THE MILLENNIUM GALAXY CATALOGUE 

The numerical resul ts presented i n thi s work use the 
galaxy catalogues of iBertone et all l|2007l ). publicly avail- 
able fr om the Millennium w ebsit j^. The Millennium simu- 
lation (|Springel et al.|[2005h uses a ACDM cosmology with 
parameters 0^ = 0.25, fib = 0.045, h = 0.73, = 0.75, 
n = 1 and erg = 0.9. The Hubble constant is parameterised 
as -ffo = 100 km s~^ Mpc~^. The Millennium cosmology 
is slightly different from that assumed by the observations, 
but the uncertainties that may derive from this discrepancy 
should be negligible in comparison with the uncertainties 
introduced, for example, by the merging time-scale. For ex- 
ample, at the highest redshifts the ratio of volumes for the 
model cosmology is 1.1% higher than the cosmology used 
in the observations. This volume effect only influences the 
calculation of the merger rates, and does not affect the cal- 
culation of the merger fractions. 

The simulation follows the evolution of about 20 mil- 
lion galaxies in a region of 500 h~^Mpc on a side. The 
galaxy catalogues were created using a variation of the semi- 
analytic model of iDe Lucia fc BlaizotI l|2007h . in which su- 
pernova feedback is imple mented using a dynamica l treat - 
ment of galactic winds (|Bertone. Stoehr fc White! !2005h . 

^ http:/ /www. mpa-garching.mpg.de/galform/virgo/millennium 
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Figure 1. The median merging time tm in units of the Hubble 
time tHubblc in the Millennium simulation as a function of the 
stellar mass of the satellite galaxy. Results are shown for z = 
1,1.5,2 and 3. 



The model improves the treatment of galaxies with stel- 
lar masses M* < 10^^ Mq , but somewhat overpredicts the 
abundance of the most massive galaxies with > 10^^ 
Mq . For a detailed description of the model, we refer the in - 
terested reader to the original paper of !Bertone et al.!l!2007h . 

The merger history of the simulated galaxies is recov- 
ered by flagging galaxies that have undergone mergers and 
by saving additional information about the time when the 
mergers occur. In order to estimate the merger fractions and 
the merger rates as closely as possible to the observations, 
we count the number of galaxies A^m that have merged at 
least once within a time-scale Tm at redshifts z = 0, 0.5, 1, 
1.5, 2 and 3. The merger fractions are then calculated as in 
Eq. ((2} and the merger rates as in Eq. (Q. 

The numerical resolution of the Millennium simulation 
allows us to correctly track the merging history of galaxies 
only for A'h > 10^'^ -^0- Accordingly, our results are not 
fully reliable below this mass threshold. Nonetheless, in a 
few cases we have chosen to display results for the merger 
rates and fractions in the interval 10^ Mq < M* < 10^'' 
Mq: the inability of the model to track the merger history 
of galaxies with 10^ Mq < < lO^ '^ Mq translates in 
to lower predicted values than in the ideal case where the 
merger history is fully accounted for. This effect is certainly 
stronger for minor merger counts, but may also affect major 
mergers. 

In the Millennium semi-analytic mode l, galaxies are 
classifled as central or satellite galaxies (jSpringel et al.! 

Central galaxies are usually the most massive galaxies 
in haloes, while satellite galaxies are smaller objects accreted 
through merging of dark matter haloes. When a smaller halo 
falls into a larger one, its galaxies become satellite galaxies 
of the resulting halo. Satellite galaxies can in turn be classi- 
fied as satellites of type 1, that is galaxies associated with a 
dark matter substructure, or of type 2, that is galaxies whose 
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dark matter substructure has completely merged with the 
central halo and can no longer be identified. 

In the model, the merging time of galaxies, tm, is the 
lifetime of a satellite of type 2 and is defined as the time- 
scale between when a galaxy loses its dark matter substruc- 
ture an d the moment when it merges with the halo central 
galaxy jSpringel et al.llioOll '). Th e merging time is given by 
the dynamical friction time-scale ijBinnev fc Tremainelll987l l 
and depends on both the virial properties of the dark mat- 
ter substructure to which the galaxy was previously attached 
and on those of the central halo. 

This definition of merging time is different from the 
"total merging time", which is instead the time between 
the moment a galaxy falls into a larger halo, and when it 
actually merges with the halo central galaxy, independently 
of when its associated substructure merges with the cluster 
halo. The total merging time is the sum of the merging time 
tm plus the time that a galaxy spends as a satellite of type 
1. Typically, the time spent by a galaxy as a satellite of type 

1 is much shorter than the time spent as a type 2. The time- 
scale for a substructure to dissolve in to a larger halo varies, 
e.g., with the parent halo and with the substructure dark 
matter mass, but on average it is shorter than 1 Gyr. As 
such, the time spent by a satellite as a type 1 is about an 
order of magnitude shorter than the merging time tm ■ 

Fig. [1] shows the median merging time of gal 
a function of stellar mass at 2: = 1, 1.5, 2 and 3. The median 
merging time tm depends only weakly on stellar mass at 

2 = 3, but at later times the more massive galaxies tend to 
have shorter merging times, on average by a factor of two, 
compared with lower mass galaxies. There is, in addition, 
a strong evolution with redshift. In particular, tm increases 
steadily with decreasing z, and it is about 3 times larger 
at « = 1 than at 2 = 3. The spikes at the highest stellar 
masses are the result of poor statistics, since only a handful 
of galaxies populate these mass bins. 

The merging time-scale in the semi-analytic model de- 
pends on the virial masses of the merging galaxies, but not 
on their baryon con tent. Therefore, because the models of 
iBertone etlo] (120071 ) and lDe Lucia fc BlaizotI l|2007i ') use the 
same merger trees, galaxies always merge at the same time 
in both models. However, since these two models do not pre- 
dict the same stellar mass, our results for the merger fraction 
and merger rates as a func tion of stellar mass a nd re dshift 
might differ from those of Lucia fc BlaizotI (|2007l ). We 
discuss the difference between the two model predictions in 
Subsection 15.11 

In the semi-analytic model the mass ratio over 
which a merger is defined as major is assumed to be 
m-sateiiite/mcentrai > 0.3, whcrc msatoiiitc is the mass in 
baryons of the merging galaxy and mccntrai is the mass of 
the central galaxy. This may not be fully consistent with 
the observations, where best estimates suggest that major 
mergers are those with mass ratios larger than ~ 0.25. This 
is quantified in Fig. (2] which shows the predicted fraction of 
mergers with mass ratios larger than 0.3, 0.25, 0.1 and 0.01 
as a function of redshift, for galaxies with > lO" M0 . 
While almost all mergers have a mass ratio larger than 0.01, 
the fraction of mergers with larger mass ratios decreases 
quickly. In particular, we find that mergers with a mass ra- 
tio larger than 0.25 are about 20 per cent more numerous 
than mergers with a mass ratio of 0.3. This implies that our 
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Figure 2. The fraction of galaxy mergers with a mass ratio larger 
than 0.3 (soUd line), 0.25 (dashed line), 0.1 (dot-dashed line) and 
0.01 (long dashed line) as a function of redshift, for simulated 
galaxies with > 10^ Mq . 

estimates of the merger fraction and merger rate may be 20 
per cent lower than the observations simply because of the 
Millennium simulation assumption of a higher mass thresh- 
old for defining major mergers. We will discuss this further 
in Sec. [5] 



4 RESULTS 

In this Section we present our results for the merger history 
of galaxies in the Millennium simulation. We first present 
results of predicted merger fractions and rates as a function 
of stellar mass in Subsection 14.11 and we then proceed to 
compare the simulation predictions with the observational 
results in Subsections 14.21 14.31 and 14.41 

4.1 Basic Features of the Predicted Merger 
History 

In this Subsection we present results on the evolution of the 
merger fraction and the merger rate as a function of stellar 
mass and redshift in the Millennium simulation. 

Fig. |3] shows the predicted evolution of the merger frac- 
tion with redshift as a function of stellar mass for both major 
(upper panels) and minor mergers (lower panels). To investi- 
gate the dependence of the merger fraction on the time-scale 
over which mergers are counted, we show results for two 
cases: one in which we have considered mergers that have 
occurred within a time-scale Tm = 0.4 Gyr (left panels) and 
a second case with a time-scale Tm ~ 1 Gyr (right panels). 
We find little difference between the results obtained using 
the two different time-scales. Both the shape and the general 
history of the minor and major mergers are very similar at 
high-z, making the number of mergers per unit time evolving 
slowly. 
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Figure 3. Predictions for the merger fractions of galaxies as a function of stellar mass and redshift in the Millennium simulation. The 
left panels show the merger rates calculated for a time-scale rm = 0.4 Gyr, while the right panels assume a time-scale rm = 1.0 Gyr. The 
merger fractions are shown separately for minor mergers (lower panels) and major mergers (upper panels). 



There is a slight difFerence between the predicted evo- 
lution of the merger fraction at high and low stellar masses, 
with a change in slope at about Af* ~ 10^°'* M0 . The ma- 
jor merger fraction is lowest at low redshift (blue solid line) 
and is highest for the most massive systems with Mi, > 10^^ 
M0 at all redshifts. At high redshifts there are few galax- 
ies with stellar masses larger than M* ~ 10^^ M0 and this 
result is affected by small number statistics. Within the ob- 
served sample, there are no galaxies with M* > 10^^ Mq 
and very few are found at > lO" '^ Mq . The predicted 
merger fraction for galaxies with M* ~ 10^'' M© at 2; = 
is /m < 0.01 for mergers within a time-scale Tm ~ 0.4 
Gyr. This result is roughly consistent with what is found 
when examining the n earby merger fraction history (e.g., 
iDe Propris et al.ll2007l . 



The minor merger fractions shown in the lower panels 
of Fig. |3] display a similar trend to the major mergers, plus 
a few additional interesting features. The first is that minor 
mergers are very common in massive galaxies, with nearly all 
galax;ies at M*, > 10^^ Mq having a minor merger occurring 
within the past 0.4 Gyr at every redshift. At lower stellar 
masses, there is an increasingly large difference in the minor 
merger fraction evolution, with a high ratio of /m, minor > 0.3 
at 2: = 3, which declines rapidly with redshift by an order 
of magnitude at 2 = 0. This means that in the Millennium 
simulation less massive galaxies undergo many more minor 
mergers at high redshift than in the recent past. 

One could argue that if the mass resolution of numerical 
simulations could be increased indefinitely, the minor merger 
fraction would be unity for all galaxies. This is certainly 
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Figure 4. The merger rate in tlie Millennium simulation as a function of stellar mass. The left panel shows the major merger rate and 
the right panel the minor merger rate. In both panels, each line shows results at different redshifts. 



true if one considers as minor mergers all accretion events, 
independently of the mass ratio between the central galaxy 
and the infalling cloud. As we show in Fig. (2] between 80 
and 90 per cent of mergers in the Millennium simulation 
have mass ratios larger than 0.01. This mass ratio is a small 
number and includes accretion of small satellites that do not 
affect the morphology of the central galaxy and would not 
be found as members of galaxy pairs. A fraction of these 
merging satellites have stellar masses below the resolution 
of the simulation, especially when galaxies with M* < 10^° 
Mq are considered. However, above M^, > 10^" Mq , most 
events we consider as minor mergers have a non-negligible 
mass ratio. Imposing a minimum value for the mass ratio for 
mergers to be considered in our analysis, as for example a 
minimum mass ratio of 0.05, would modify the minor merger 
results by about a factor of 2 for the most massive galaxies, 
but would not change the qualitative behaviour of the minor 
merger fractions and rates. The redshift evolution of the 
minor merger fraction as a function of stellar mass shown in 
Fig. [3] is an indication that resolution effects alone do not 
shape the distribution of the minor merger fractions as a 
function of stellar mass. According to Fig.O major mergers 
represent between 5 and 10 per cent of all merger events in 
the Millennium. 

Fig.|4]shows predictions of the merger rate as a function 
of stellar mass in units of Mpc~^ Myr~^. The evolution of 
the merger rate with redshift is shown for both major merg- 
ers (left panel) and minor mergers (right panel). The differ- 
ent lines indicate results at different redshifts. The merger 
rate tends to be higher for less massive galaxies and lower 
for the most massive galaxies. The shape of the merger rate 
as a function of stellar mass mostly reflects the shape of the 
stellar mass function ngm(z), which enters the definition of 
the merger rate according to Eq.|4l In fact, while the number 
density of galaxies varies by up to five orders of magnitude 
for lO'' M© < M* < 10^^ Mq systems, the merger fractions 



vary by at most a factor of 100, as shown in Fig. [S] The 
major and the minor merger rates are fairly constant for 
galaxies with Ah < 10^^ Mq , at 2 > 1 and decrease with 
time for z < 1. The merger rates of galaxies with Ah > 10^^ 
Mq strongly decrease with redshift. This partially reflects 
the increase in the number density of galaxies with time. 



4.2 Comparing the Observed and Predicted 
Merger Fraction Evolution 

In this Subsection, we examine in detail how the observed 
merger fractions compare with those predicted by the Mil- 
lennium simulation. The merger fraction, as defined in Sub- 
section 12.11 is the fraction of galaxies undergoing a merger 
within a given stellar mass and redshift range. 

Fig. [5] shows a comparison between the observed and 
simulated merger fractions for galaxies with > lO^'^ Mq 
(upper panels) , M^, > 10^° Mq (middle panels) and 10^ Mq 
< Ali, < 10^" M0 (lower panels) as a function of redshift. 
To demonstrate the effect of the time-scale for merging, the 
simulation results in the left panels are for a time-scale Tm = 
0.4 Gyr, those in the right panels for Tm = 1 Gyr. 

The predicted major merger fractions for the most mas- 
sive systems with M* > lO^'^ Mq , shown in the upper pan- 
els of Fig. O display a good agreement between observations 
and simulations at z < 2. There is, however, a large differ- 
ence between the models and the data at z > 2, where the 
observed merger fraction is about 5 to 10 times higher than 
predicted by the model. 

The situation is rather different for systems with Af* > 
10^" Mq , for which the predicted major merger fractions are 
systematically lower than the observed values by a relatively 
large factor, independently of the time-scale considered. The 
predicted major merger fraction for a time-scale Tm — 0.4 
Gyr is about an order of magnitude smaller than observed 
at 2 ~ and the difference increases significantly at higher 
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Figure 5. Comparison between the observed and simulated merger fractions for galaxies with Af* > 10^^ M0 (upper panels), Af* > 10^'^ 
Mq (middle panels) and 10^ Mq < Af* < 10^'' Mq (lower panels). The dashed and solid lines indicate the minor and major merger 
fractions respectively. Predictions are shown in the left panels for a time-scale of Tm = 0-4 Gyr and in the right panels for Tm = 1 Gyr. 
The observed data points, described in Section |2] are filled circles for CAS results and diamonds for pair fraction results. 
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Figure 6. Comparison between the observed values of tlie inverse of the merger rate per galaxy F and predictions from the Millennium 
simulation as a function of redshift. The value of F is roughly the time-scale between mergers. The predicted value of F is calculated 
assuming Tm = 0.4 Gyr. Results are shown for two mass bins: Af* > 10^'^ Mq (left panel) and Af* > 10^^ Mq (right panel). 



redshifts. It is most likely a coincidence that in this mass 
range the observed merger fractions nearly match the pre- 
dicted minor merger fractions. The minor merger fraction 
changes slope at 2: ~ 1, with a fast decline with decreasing 
redshift at 2: < 1 and a slow increase with redshift at z > 1. 

As already mentioned, the largest uncertainty in esti- 
mating merger fractions in simulations is the uncertainty in 
the time-scale Tm to which the observations are sensitive. 
Fig. [5] shows that the observed major merger fractions are 
a factor of a few higher than in the Millennium even for a 
time-scale of Tm = 1.0 Gyr. This is an indication that in 
this mass range either there are not enough major mergers 
in the simulations, or a larger number of minor mergers are 
counted as major in the observations. Alternatively, if the 
time-scales we assume are incorrect, they would have to be 
Tm = 4 — 10 Gyr to match the predictions. 

Finally, in the bottom panels of Fig. [5] we show the 
merger fraction evolution for galaxies with 10^ Mq < Af* < 
10'" Mq . As we mentioned earlier, the simulated galaxies 
with stellar masses lower than about 10^'^ Mq are affected 
by resolution and their merger history might not be fully 
resolved. Consequently, the model results shown here most 
likely underestimate the real merger fractions. With this in 
mind, the observed major merger fractions are about an or- 
der of magnitude higher than the predicted ones, confirming 
what is seen for the mass range Af* > 10^° Mq . However, it 
is worth pointing out that the overall merger fraction evolu- 
tion for systems in both mass ranges has a similar shape as 
the observed one. We discuss this further in the next Section. 

The merger fractions at z > 2, that is the two data 
poin ts at the h ighest redshifts in all ranges, are taken 
from lBluck et al.1 12009. ), who measure merger fractions from 
galaxy pairs in the GOODS NICMOS Survey. The estimate 
of the merger fraction with galaxy pairs is somewhat less 
likely to be biased by misinterpreting the nature of the 
merger than with morphological methods. However, other 



uncertainti es might affect the estimate of merger fractions 
usmg pairs. iKitzbichler fc White! l|2008l ) argue that the time- 
scale for merging generally considered by observers is under- 
estimated by at least a factor of two. This means that the 
two high redshift data points might be overestimating the 
merger fractions by a similar factor. Since the difference be- 
tween the observed and predicted merger fractions at z = 2.5 
is about a factor of ten, it is likely that there is an additional 
discrepancy that remains unexplained even after taking into 
account the correction suggested by IKitzbichler fc White! 
(,2008. ). Changing the merging time-scale from 0.4 Gyr to 
1.0 Gyr, as in Fig. [5] does not help solve the discrepancy, 
as the predicted merger fractions are then too low at high 
redshift and slightly too high at low redshifts. 

A few studies l|Conselice et al.ll2008l : Isiuck et al.ll2009l ) 
have recently investigated the behaviour of the inverse of the 
merger rate per galaxy, defined as F = Tm//gm. The merger 
rate per galaxy is equivalent to the merger rate without the 
dependence on the stellar mass function and is basically the 
average number of mergers a galaxy experiences per unit 
time. As such, F represents the average amount of time be- 
tween mergers, and is independent of the merging time-scale 
Tm in the models. We show results for F in Fig.|6]for the stel- 
lar mass ranges A//* > 10'" Mq and Ah > lO" Mq . 

In all the stellar mass bins we consider in the models, F 
is high at high redshift, declines to a minimum at around z ~ 
1.5 and then rises again at lower redshifts for major mergers. 
In other words, the predicted time between mergers at 2 > 2 
is longer than it is at z ~ 1.5. The predicted F steadily 
increases with decreasing redshift at z < 1.5, for all mass 
ranges and for both major and minor mergers. On the other 
hand, at 2 > 1.5, the predicted F slowly increases for major 
mergers, and decreases for minor mergers. T his is consistent 
with observations of massive galaxies (e.g., [Conselice et al.l 
I2OO8! ). The predicted F for major mergers is about ten times 
larger than observed for galaxies with Af* > 10'° Mq , while 
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it matches the observations well at low redshifts for galaxies 
with Mi, > 10^^ Mq at z < 2. This is a further indication 
that the merging history of galaxies with Af* < 10^^ do 
not agree in simulations and observations. 



4.3 Simulated Merger Fraction Parameterisations 

A common thing to try when examining the evolution of 
the merger fraction is to pa rameterise the evoluti on in 
various ways. As discussed in IConselice et al.l (|2008l ) . the 
two most popular ways to parameterise the merger frac- 
tion are through a power-law and through a combined 
exponential/power-law Press- Schechter-like function. The 
power-law fit is by far the most common method, but is 
becoming less fashionable as it appears th at the observed 
merg er fraction turns over at high redshifts (|Conselice et al.l 
|2008| ). while a power-law continues to increase at all red- 
shifts. The form of the power-law evolution is given by: 



U{z) = /o X (1 + zY 



(5) 



where /m(z) is the merger fraction at a given redshift, /o is 
the merger fraction at z = and m is the power-law index 
for characterising the merger fraction evolution. Investigat- 
ing how to parameterise the increase in the predicted merger 
fraction evolution is another way to determine how the ob- 
served and predicted merger fractions differ. We find from 
previous work that the index m on this power-law increase 
in the merger fraction is typically m ~ 2 — 3, which has been 
found using a variety of t echniques (e.g., IConselice et al.l 
I2OO8I : IConselice et al.ll2009l and references therein) . 

We find that the merger fraction in the Millennium sim- 
ulation can be characterised by a power-law index m — 1.6 
for 2^1 and m = 0.99atl<2;^3 for galaxies with stellar 
masses > 10^^ Mq , and m = 3 for galaxies with masses 
M* > 10^" M0 at z < 1. These power law fits are only valid 
up to the merger fraction turnover at 2 ~ 1 — 1.5. While the 
amplitude of the observed merger fractions are higher than 
the simulation, the slope is fairly similar to the predictions. 

Another way of characterising the merger fraction evo- 
lution, which is based on th eoretical argum ents that use the 
Press-Schechter formalism (|Carlberg| Il990l ). is a combined 
power-law exponential evolution of the form: 



U = a{l + zr xe^p[p{l + zf] 



(6) 



This combined power-law/exponential description repro- 
duces the obser vations better than a simple power-law 
l|Conselicell2006al ) and fits the merger fractio n predictions, as 
well as the data (e.g.. IConselice et al.ll2009l ). The behaviour 
of the merger fraction evolution can be interpreted as either 
the result of mergers occurring later for massive galaxies 
in haloes, due to the dynamical friction time-scales, or to 
the fact that there are not many very massive galaxies to 
merge with at early times, resulting in a lowered merger 
fraction. In general, it appears that all galaxies, with the 
exception of the most massive ones with > 10^^ Mq , 
have a turnover in their merger fraction history and can be 
fit by an exponential/power-law, which is likely the correct 
form for parameterising the merger fraction history. 



4.4 Merger Rate Evolution 

In this Subsection we present results for the merger rate 
in the simulations and compare them to the observational 
data. 

The merger rate ^(2) is currently the most uncertain 
merger quantity to measure in observations. Estimates will 
likely become more straight-forward as our understanding 
and knowledge of galaxy number densities, merger fractions 
and merger time-scales improve. On the other hand, the 
merger rate is a relatively easy quantity to measure in sim- 
ulations, once a proper definition has been agreed upon. 
Since ^{z) is calculated by dividing the galaxy merger frac- 
tion, which scales with the time-scale, by the time-scale, the 
merger rate is independent of the time-scale Tm used to mea- 
sure the merger fractions. This is true under the assumption 
that the merger rate slowly varies with time and can be con- 
sidered approximately constant within a short time-scale. 
However, the assumption breaks down for large time-scales, 
because the merger rate itself is not intrinsically constant 
over large time spans. 

The comparison of the predicted and observed merger 
rates is shown in Fig. [7] The predicted merger rates for the 
10^ M0 < M* < 10^° Mq galaxies (upper left panel) agree 
well with the data, despite the shortcomings due to lack 
of resolution in the simulation. The results for the other 
mass bins, that is M* > 10^" Mq (upper right panel) and 
M* > lO^'^ Mq (lower panel), agree with the observations 
within the error bars. 

The merger rates display a surprisingly good agreement 
between the observations and the Millennium simulation, 
which is not seen for the merger fractions of galaxies with 
Af* < lO^'^ Mq . Since the observed and predicted merger 
fractions agree only for the most massive galaxies, the agree- 
ment for the lower mass bins is likely coincidental. As dis- 
cussed in Subsection l2.3l the merger rate is a function of the 
galaxy merger fraction /gm, the number density of galaxies 
rigni and the time-scale for merging Tm. Any one of these 
quantities can affect the measurement of the predicted and 
observed merger rates. The merging time-scale is constant 
in our model, but we cannot exclude an inconsistency with 
the time-scale assumed by the CAS and pair methods. How- 
ever, this reason alone is unlikely to explain the discrepancy 
between the merger fractions and rates in the lowest mass 
bins, given the agreement at the highest masses. 

The one quantity we have not examined yet is the pre- 
dicted number density of galaxies within the Millennium 
simulation. Since we are selecting galaxies through a stel- 
lar mass cut, the predicted number density of these sys- 
tems affects the final values of the merger rates. We examine 
this in Fig. [S] which shows a comparison between the Mil- 
lennium predictions and observations for the number den- 
sity of galaxies with Af* > 10^^° Mq (filled triangles) and 
Mj, > 10"*^"*^ Mq ( empt y squares) . Data points are taken from 
IConselice et al.l (120071). The so lid line shows results for the 
model of Bertone et al.l |2007^ ■ the dashed line for that of 
iDe Lucia &: Bl aizot' (200 7|). 

The model of Bert one et al. roughly predicts the 

number density of galaxies w i th Mj , > 10^^ Mq , while the 
model of iDe Lucia fc BlaizotI l|2007l ) underpredicts it at all 
redshifts. Both models overestimate the observed values for 
galaxies with Af* > 10^° Mq by almost an order of magni- 
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Figure 7. Comparison between the observed merger rates in units of co-moving Mpc per Myr and predictions from the Millennium 
simulation. The different lines are the same as in previous figures and show the major merger rate (solid line) and the minor merger rate 
(dashed line) in three different mass ranges: 10® Mq < < 10^'' Mq (upper left panel), > 10^'' Mq (upper right panel) and 
Af* > lO^l M0 (lower panel). 



tude. This results in the high value derived for the merger 
rate of galaxies with Af* < 10^^ M0 , since, as we know 
from Fig. O the predicted merger fraction is low compared 
to the data. We therefore conclude that the apparent agree- 
ment seen in Fig. [7] is partly a coincidence for the stellar 
mass interval M* > 10^'' M0 . We believe that the agree- 
ment is instead genuine for the more massive galaxies with 
> 10" Mq . 



5 DISCUSSION 

Tracing the galaxy merger history observationally is just 
now becoming a mature enough area of astronomy that 



it can be compared with theory. The modern concept for 
how galaxies and structures form in the universe is based 
on the idea that galaxies accrete a significant fraction of 
their mass through mergers. While the ACDM paradigm 
predicts many features of the observed galaxy population 
successfully, many galaxy properties are difficult to explain 
and reproduce in models. These include the mass profiles o f 
galaxies, the missing satellite problem (jMoore et al.|[l999h 
and the fact that at high re dshift there are more massive 
galaxies than predicted (e.g., IConselice et al ] l2007l ). These 
differences might be due to the way that the physics of star 
formation is implemented in the models, or to the underly- 
ing dark matter halo history, or perhaps even the underly- 
ing cosmology. Comparing the merger histories predicted by 
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Figure 8. The observed stellar mass density as a function of red- 
shift for galaxies with M* > 10^° Mq (filled triangles) and Ah > 
10^^ M p) (empty squares). Data points are fromlConsclicc ct al. 

2007). The solid line shows results for the mod el of Bcrtonc ct al. 

20071 ). the dashed line for that of iDe Luciafc Blaizot (,2007.) . 



models of galaxy formation with available data is a new way 
to help test our understanding of galaxy formation. 

As discussed in Subsection 14.41 we find that major 
merger rate predictions match the observations fairly well. 
The merger fraction predictions discussed in Subsection 14.21 
match the observations well for galaxies with AI-^ > 10^^ 
Mq at 2 < 2, but diverge significantly when less mas- 
sive galaxies are considered. For these galaxies, the ob- 
served major merger history is more substantial than pre - 
dict ed by the semi-analytic rn odels of lBertone et al.l (|2007l ) 
and lDe Lucia fc Blaizot |20o3). Previous investigations have 
highlighted a similar problem when examining galaxies in 
pairs compared with the Millennium simulation ( Mateus 
l2008l ). 

In the following we discuss several possible explanations 
for the discrepancy between the observed and the predicted 
merger rates and fractions of galaxies with intermediate 
masses, including observational biases in finding mergers, 
problems with the underlying mass assembly history and 
the predicted star formation history of galaxies. 



5.1 Variations due to the physical modeHing of 
star formation and feedback 

One major issue we have not addressed yet is the fact that we 
are selecting galaxies by their stellar mass. The merger frac- 
tion and rate vary as a fu nction of stellar mass, both within 
observed galaxies (e.g., IConselic^ l2003al : IConselice et al.l 
[2OO8. ) and within models (see Fig. |3]). An uncertainty in the 
determination of the stellar mass of galaxies could translate 
in to an uncertainty in the measured merger values. Uncer- 
tainties in the star formation history would place galaxies 
into too high or too low stellar mass bins and affect the mass 



ratio of the merger event. This in turn afi'ects the predicted 
major merger fractions and rates. 

In the semi-analytic model, galaxies that fall into a 
larger halo cannot form many stars after they become satel- 
lites, as they cannot accrete gas from their environment 
and quickly consume or expel the interstellar medium they 
might have possessed when they were central galaxies them- 
selves. This is a well-known problem within semi-analytic 
models and is one of the causes of the excess red popula- 
tion of low mass galaxies on the col our-magnitude diag ram, 
which are not seen in observations l|Bertone et al.ll2007l ). As 
a consequence, a larger fraction of predicted merging galax- 
ies might have lower mass ratios than what is necessary to 
be counted as major mergers, and a larger number of major 
mergers could end up being counted as minor mergers. A fac- 
tor that could counteract this effect, though, i s stellar mass 
loss through tidal strippin g in cluster galaxies (Bullo ck et al.l 
l200ll : [Murante et al.ll2007l ). which is currently not accounted 
for in the semi-analytic model. 

As ment ioned in ^ the mer ger fractions and rates in 
the model of lBertone et al.l (j200'/) may differ from those of 
other models based on th e Millennium simulation, in cluding 
for example the model of iDe Lucia fc Blaizoj l|2007l ). These 
two models differ only in the treatment of supernova feed- 
back, but otherwise use the same merger trees and the same 
treatment of gas physics. Star formation, AGN feedback, 
gas cooling, merging and other relevant processes are done 
in the same way, using the same parameters. 

The SN feedback scheme is different in the following 
ways. In the|De Lucia & Blaizot (2007,) model the SN feed- 
back is based on simple formulae of e nergetic balanc e , as in 
other semi-analytical models. In the iBertone et al.l (|2007l ) 
model, the evolution of galactic winds is followed by solving 
the equations of motion of the outflows, which are mod- 
elled as pressure-driven cosmological blastwaves when they 
emerge from galaxies and become momentum-driven snow- 
ploug hs when cooling sets in . The treatment of SN feedback 
in the iBertone et al.l ||2007|) model implemen t s a fa ster re- 
cycling of gas than in the loe Lucia fc BlaizotI (|2007l ) model, 
which suppresses star formation more efficiently in dwarf 
galaxies and less efficientl y in massive galaxies than the feed- 
back prescriptions in the lOe Lucia fc BlaizotI (|2007h model. 
This results in a larger number of massive galaxies and 
in a shallower s lope at the fa i nt end of the stellar mass 
function in the Bertone et al.l (|2007l ) model than in the 
iDe Lucia fc BlaizotI (|2007l ) model. 

In Fig. [9] we quantify the difference between these 
two models by showing the ratio of the predicted major 
merger rates ( l eft pa nel) and fractions (right p anel) of the 
IBertone et"al] (|2007l ) and iDe Lucia fc BlaizotI (|2007l ') mod- 
els. Results are shown as a function of stellar mass for the 
redshifts 2 = 0, 1,2 and 3. The dashed horizontal line indi- 
cates where the model predictions are equal. The two mod- 
elled merger rate predictions differ by no more than a factor 
of three, in the most extreme case, with the largest difference 
being in the values of the merger rates at 2 > for galaxies 
with M* < 10^^ M0 . As discussed above, a possible expla- 
nation for this difference is the h igher abundance of galaxies 
in this stellar mass range in the lDe Lucia fc BlaizotJ (|2007| ) 
model. The predicted merger fractions are roughly similar, 
with the important exception o f a narrow stellar m a ss bin 
around A'L ~ lO" Mq . In the iDe Lucia fc BlaizotI (|2007l ) 
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Figure 9. Comparis on between the major merger rates (left panel) and major merger fractions (right panel) predicted by the BDT07 
l|Bertone et al.ll2007^ and by the DLB07 l|De Lucia fc BlaizotI EoOTI ) models as a function of galaxy stellar mass. Results are shown 
for redshifts z = 0,1,2 and 3. In both panels, the y— axis shows the ratio of the values predicted by th e two models. The dashed 
horizontal line indicates where the model predictions are equal. Both figures demonstrate that in general the lPe Lucia &: Blaizo 3 l|2007l ) 
model predict s high er values for the merger rates and fractions for M* < lO^'^ — lO'^^'^ Mq at almost all redshifts. Conversely, the 
iBertone et al.l l|2007l ) model provides higher merger rates and fractions in the highest stellar mass bins. 



model the merger fractions for the Af* > 10 Mq galaxies 
nearly match the observations if the time-scale for merging 
is Tm = 1 Gyr. Assuming Tm = 0.4 produces results that 
slightly underp redict the data. The predic ted major merger 



slightly underp redict the data. 1 he predic ted major merger 
fractions in the lDe Lucia fc BlaizotI (12007 ) model differ by a 
very minor amount from the Bertone et al.l l|2007h model for 
galaxies at M* < 10^" Mq . As a consequence, the different 
feedback prescriptions used in the two models do not signif- 
icantly help to better match the observations at low stellar 
masses, where the discrepancy is largest. 



5.2 Variations in the Merger Histories 

The discrepancy between the observed and predicted major 
merger fractions could also result from a lack of merging 
events in the simulation itself, or to systematic differences 
in the calculations of merger fractions, or in the definition 
of what is a merger. 

One difference between the model and the observations 
that we are aware of, and which we discussed in Section [S] 
is the definition of what is a major merger. Observations 
define mergers as major when the mass ratio between the 
merging galaxy and the central galaxy is greater than 1:4. 
In the Millennium simulation, major mergers happen when 
the mass ratio is greater than 0.3. This difference may ac- 
count for some of the discrepancy in the estimates of the 
merger fractions and rates: if a mass ratio 0.25 were used 
also for the simulations, the estimate of the merger fractions 
would increase by about 20 per cent, as it would also include 
mergers that are currently not counted as major. This would 
improve somewhat the agreement between the observed and 
predicted merger fractions in some stellar mass bins, but it 



would not significantly improve the overall agreement and 
it would not solve the discrepancy at z > 2. iMateusI (|2008| ) 
also points out that the observed merger rates and fractions 
are better reproduced by the models when mergers with a 
mass ratio smaller than 1: 4 are inc l uded in the estimation 
of the merger rate. Indeed, iMateusI l|2008l ) best match with 
the observations is achieved for a mass ratio of 1:10. 

It is also not entirely clear what is the CAS major 
merger ratio sensitivity. This is generally assumed to be 0.25, 
but it could equally be 0.3, or lower. We know this is unlikely, 
base d on c omparin g CAS and pairs for nearby galaxy merg - 
ers (|Hernan dez-T oledo et alllsioOSl : [Pe Propris et aLll2007l 'i. 
For a basic agreement at high redshift or at low stellar 
masses, the CAS method would have to be sensitive to merg- 
ers with mass ratios down to 1:10. This is again unlikely to 
be the major reason for the difference between models and 
observations, because the variations involved are too small 
in comparison with the discrepancies in the results. How- 
ever, it is a further indication that the mass ratio threshold 
to identify major mergers is indeed a rather volatile factor 
both in the observations and in the simulations. 

Finally, the mismatch between the observations and 
simulation results might be due to the inability of the Millen- 
nium simulation to produce enough major mergers. This can 
happen, for example, if there are too few galaxies in haloes, 
or if the predicted merging time-scale is too long. The lack 
of galaxies in haloes can be tested by measuring the bias of 
massive galaxies in the Millennium simu lation, which have to 
be more clustered than the dark matter. [Zehavi et al.l (|2005l ) 
found that the predictions of semi-analytical models agree 
fairly wel l with the cluster ing properties of nearby galaxies. 
However, ICoil et al] l|2008l ) found that at high redshift sim- 
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ulated blue galaxies are much less clustered than real blue 
galaxies, especially at small scales. As the vast majority o f 
merging galaxies at 2: > 1 are blue (jConselice et al.ll2003lJ ). 
this is a potential problem of the models, which translates in 
to an underestimate of the merger history. In fact, a lack of 
a strong bias, or equivalently a lower halo occupation num- 
ber, implies that massive systems do not experience as many 
merger events as they would if they were more clustered. 



5.3 Comparison with Previous Work 

In this Subsection we briefly compare our findings with re- 
sult s from previou s work s. 

IGuo fc Whitel (|2008l ) examine how galaxies grow within 
the Millennium model. They find that the most massive 
galaxies have a formation history dominated by major merg- 
ers, while f or less massive galax i es the same is true only 
at 2 > 1. iKitzbichler fc Whitl (|2008l ) use the model of 
iDe Lucia &: BlaizotI ( 2007 ) to calibrate the relationship be- 
tween the abundance of close galaxy pairs and the rate 
of galaxy mergers at high redshifts. They flnd that close 
galaxy pairs merge within a few Gyr, and that this merging 
time-scale only weak l y inc reases with decreasing redshift. 
IKitzbichler fc Whitel (I2OO8I ) also argue that while the use 
of close pairs for estimating the merger history of galax- 
ies is a reliable tool, the merging times used to calculate 
the merger rates in observations are at least a factor of two 
too short with respect to their findings, which translates in 
to an overestimate of t he merger rates by a similar factor. 
iPatton fc Atfieldl (|2008l ) show that this is not the case when 
an appropriate correction factor is included in the equa tions. 
Similar conclusions are reached bv lWetzel. eTail l|2008l ). who 
find that only pairs at very small separations can be consid- 
ered a reliable proxy for the glo bal merger pop ulation. 

Using galaxy pair fractions, iMateud l|2008l ) finds a sim- 
ilar qualit ative behaviour with redshift for the merger rate 
as we do. iMateud (I2OO8I) also compa res the results from 
the model of'Pe Lucia fc Bl aizotI (120071 ') with those from the 
model of iBower et a l. (2006i) and finds that the lBower et al.l 
(|2006l ') model better reproduces the qua litative behaviour oi 
the m erger rate with redshift than the iDe Lucia fc BlaizotI 
l|2007l ) model, although the predicted values a re still lower 
than observed. This is interesting, because Bo wer et al 



200£) use different merger tre es than iDe Lucia fc Blaizot 
20071 ) and lBertone et al.1 (|2007l ) and have different prescrip- 



tions for merging that might better reproduce the evolution 
of the observed galaxy population. 

Other studies have compared the obse rved galaxy 
merger properties with the Millennium model JCenel et al.l 
l2008l : iPatton fc Atfieldll2008l ) . iPatton fc Atfieldl (|2008l ) com- 
pare the number of close galaxy pairs (with separations 
smaller than 20 h^^ kpc) in the simulation and in a sample 
of SPSS galaxies at low redshift and find good agreement. 
IPatton fc Atfieldl (|2008l ) estimate that at least 90 per cent 
of major m e rgers occur between galaxies fainter than L*. 
iGenel et al] ^00§) flnd that there are not enough major 
mergers in the Millennium simulation to transform the pop- 
ulation of z = 2 galaxies in to the elliptical galaxies observed 
at z = 0, which implies that in the model internal evolution 
must play a more important role to shape their properties 
than in the real universe. 



6 CONCLUSIONS 

We have carried out a comparison between galaxy merger 
history predictions based on the Millennium simulation and 
a collectio n of observed galaxy merger fra ctions and rate s 
taken from Conselice et al. ^ (2003b), Conaelice et all (|2008l ). 



iBluck et al 



(I2OO9I ) and IConsehce et all ()2009V The ob- 



served merger rates and fractions are measured using the 
CAS system and by measuring the number of galaxy pairs 
at 2: > 2. 

Our main results for the Millennium galaxy formation 
model of lBertone et all (|2007l ) can be summarised as follows; 

• The predicted major merger rates and fractions vary 
signiflcantly with time and stellar mass. Massive galaxies 
experience on average more frequent major merger events 
than less massive galaxies, although the merger rate in gen- 
eral is higher for less massive galaxies. 

• The model predicts that the major merger fraction in- 
creases with redshift between z — and about z ~ 1.5, 
reaches a peak at about z ~ 1.5 and then declines at higher 
redshifts. This result is independent of the time-scale Tm 
over which mergers are identified and is common to all the 
stellar mass ranges we consider. This is roughly consistent 
with the behaviour of the observed data at z < 2, but not at 
higher redshift, where the observed major merger fractions 
continue to increase or remain high at larger redshifts. The 
predicted values for the major merger fractions match the 
observed ones for galaxies with > 10^^ Mq , but under- 
estimate them when less massive galaxies are considered. 

• The major merger rates predicted by the Millennium 
simulation roughly agree with the observed rates in all stel- 
lar mass ranges and reproduces well their evolution with 
redshift. However, since the major merger fractions agree 
with observations only for galaxies with > 10^^ M© , 
the agreement at the lower stellar masses is likely coinciden- 
tal. 

• We discuss several issues that might introduce uncer- 
tainties in the measurement of the merger fractions in the 
models and in the data. Some disagreement can be alleviated 
at some redshifts and stellar mass ranges, but becomes worse 
in other cases. Among other things, a different implementa- 
tion of SN feedb ack in the semi-analytic mo del, such as that 
implemented bv lDe Lucia fc BlaizotI (|2007| ). produces varia- 
tions in the results by up to a factor of a few. The mass ratio 
to which different techniques are sensitive might also be a 
source of uncertainty, but probably only a minor one. Like- 
wise, the stellar masses in the Millennium simulation may be 
affected by processes, such as strangulation and tidal stellar 
stripping, that are not realistically implemented, and that 
ultimately affect the mass ratio and the number density of 
mergers. It is also possible that models are not producing 
enough major mergers in galaxies with M* < lO^'^ Mq . 

The fact that the merger predictions do not match the 
observed major merger history is potentially an explana- 
tion for other known problems in the semi-analytic models 
associated with the Millennium simulation. These include, 
for example, the overproduction of red, low mass galaxies in 
colour-magnitude diagrams. This problem might be partially 
alleviated if the galaxy merger history is more substantial 
than predicted. 
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